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I. INTRODUCTION 



The UV complete theory of gravity has been a long-felt want since general relativity (GR) 
and quantum field theory were established. For past two decades, string theory has been 
considered to be a strongest and possibly unique candidate. Recently Hofava proposed a 
more practical/economical way to achieve this goal based on quantum field theory by giving 
up the general covariance of GR [H, H|. In this theory, the time and the space show different 
scaling behaviors at the UV regime. It modifies the high momentum behavior of the graviton 
propagator and renders the theory power-counting renormalizable. 

This novel feature of Hofava-Lifshitz (HL) gravity attracted much interest and various 
aspects of the theory have been investigated, including the properties and the consistency 
of the theory 0, 0, @, 0], the properties of classical solutions including black holes and 
their thermodynamic properties [§1, 0|, the cosmological aspects (lol . the perturbation 
spectrum and the gravitational wave production [lij, and the phenomenological sides (l3| . 
One notable feature of the theory is that the general covariance (full 4D diffeomorphism) 
is reduced to the foliation-preserving diffeomorphism. Thus, we have additional degree of 
freedom, a dynamical scalar 0,0], but with a first-order (in time) equation of motion. It gives 
rise to subtleties in canonical quantization, instability problem, and strong coupling problem 
(ij. On the other hand, the theory must be reduced to GR in the IR limit. It means that 
the symmetry must be enlarged at IR, the so-called emergent symmetry. Since HL gravity is 
supposed to be a quantum theory, it depends on the running of coupling constants, of which 
we do not have full understanding yet. This may lead to phenomenological difficulties, for 
example, the energy dependence in the limiting speed c 2 (E) and Sc 2 (E). To get the better 
control over these problems, we may impose the projectability condition on the metric and 
the detailed balance conditions on the coupling constants of the theory. So far, it seems that 
imposing the detailed balance condition up to the IR regime is problematic and imposing 
the projectability is favorable. 

Though the explicit quantization and the proof of the renormalizability are yet to be 
explored, it is interesting to see the consequences of such an approach at the IR regime, 
that is, in the classical solutions and to find the observable signatures through them. In this 
paper, we obtain the static spherically symmetric vacuum solutions of HL gravity imposing 
the detailed balance conditions only in the UV limit. Thus specific relations between pa- 
rameters in the four or less spatial derivative terms are not imposed. For HL gravity, the 
choice of coordinates is important because the theory lacks the general covariance and is 
intrinsically non-relativistic. To define a theory, we must specify a frame where the action 
is defined. Without knowledge about the preferred frame a priori, we choose two coordinate 
systems, the Painleve-Gullstrand (PG) coordinates and the Poincare coordinates for com- 
parison. The former coordinate system is better motivated for HL gravity since it satisfies 
the projectability condition The comparison of the solutions in two coordinate system 
will contrast the differences between projectable and non-projectable HL gravity and shed 
some light on the consequences of imposing the projectability condition. Specifically, for 
the obtained solutions without matter distribution in this paper, the former shows a long 
range effect like a surplus or deficit solid angle irrespective of the quartic derivative terms 
but the latter only has a short distance correction like the change of event horizons. In the 
presence of quartic derivative terms, the effect of those in the former is comparable to the 
electrostatic field of a point charge but, in some parameter range, the square of it has a 
negative value. 
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This paper is organized as follows. In Section [TT] we review HL gravity, setting up our 
notations. The static spherically symmetric vacuum solutions are obtained in the PG coor- 
dinates in Section HTT| and in the Poincare coordinates in Section HVl We conclude in Section 

El 



II. HORAVA-LIFSHITZ GRAVITY 

The HL gravity has the invariance under the foliation-preserving diffeomorphism t — > t(t), 
x l y x*(t, x), and a scaling behavior in the UV limit t — > £ 3 1, x l — > ix % . The action for the 
HL gravity is best described using the following dynamical variables: the lapse function N, 
the shift functions JVj, and the three-dimensional spatial metric gfy, with which the metric 
takes the ADM form 

ds 2 = -N 2 dt 2 + gij (dx i + Wdt) (dx j + N j dt) , (1) 

where N % = g^Nj. 

In the UV regime there are 5 independent sixth order spatial derivative terms saturating 
the z = 3 anisotropic scaling [Ej, 

R 3 , ROR, RijR^ R, RijOR^^ , R^jR* ^R?^ , (2) 

which make the theory less predictable. Here Rij is the Ricci tensor of g^, R = g^Ry, and 
square(D) denotes the Laplacian with respect to the spatial metric g^. Since the detailed 
balance condition spoils the IR dynamics of HL gravity and the physical motivation for its 
introduction is not manifest yet, we impose it only in the vicinity of UV fixed point. Then 
the UV dynamics in the limit of z = 3 scaling is assumed to be governed by the quadratic 
time derivative terms and square of the Cotton tensor [l|, 

ikl jm / i \ 

C lJ = — y — [R lm - -g lm Rj , (3) 

where e* jfc is the antisymmetric tensor density with e 123 = 1 and semicolon(;) denotes the 
spatial covariant derivative. Then the HL gravity action under consideration is given by 



S UL = J dtd 3 xNy/g 



a (KijK ij - XK 2 ) +£R + <j 

(4) 



pijk „lm 

+ PC^CV + 1 —^R il R km . j + c, /»',,/.'" + r]R 2 
V 9 



where Kij is the extrinsic curvature 



2N V dt 

The action (j4j) possesses 8 parameters, a, X, £, a, (3, 7, £, 77, and, in the IR limit, a, ^, 
and a terms dominate over higher derivative terms. To recover GR, the renormalization 
group (RG) running toward the IR limit must lead to 

1 c cA 

A = 1 ' a= 16^' a = ~^ (6) 

where c is the speed of light, G is the Newton's constant, and A is the cosmological constant. 
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III. STATIC SPHERICALLY SYMMETRIC SOLUTION IN THE 
PAINLEVE-GULLSTRAND COORDINATES 



The projectability condition states that the lapse function is a function of the time only, 
that is, N = N(t) in the metric (JT|). Then the time reparametrization, a symmetry transfor- 
mation in HL gravity, always allows the fixation of N = 1. The difficulties without imposing 
the projectability condition have already been discussed for quantization of HL gravity jzj, 
and then we study the so-called projectable HL gravity in this section. This projectable 
version is a different theory from non-projectable HL gravity since the foliation-preserving 
diffeomorphism cannot turn one into the other. 

To get the static spherically symmetric vacuum solution under the projetability condition, 
let us consider the static spherically symmetric metric in the PG coordinates 



ds 2 



-dt 2 



\dr + n(r)dt} 2 + r 2 (d6 2 + sin 2 9d<j) 2 ), 



(7) 



which satisfies the projectability condition. The lapse function and the shift function for the 
metric ([ID are unity and n/f, respectively. 

The action (jl]) in terms of / and n is given by 



S'hl = 47r / dt dr —j= < 



(A - 1) (2n> f 
f 



4 



n 



2(C + 277) 



2A / 


'2n' 


r \ 


, n 


-I) 2 





/ 



(rf + f-l) + a , (8) 



and the action for the matter-field is 



f r 2 
S M = 4ir J dtdr — £ M (n, /,$), 



(9) 



where $ stands for all the matter fields of our consideration. Since all the components of 
Cotton tensor vanish, C* J = under this PG metric (ED), the contribution from the sixth and 
fifth order spatial derivative terms disappears in the action and so will do in the equations 
of motion. It means that the spherically symmetric solutions of our interest are appropriate 
in describing up to the intermediate energy scale involving quartic spatial derivatives and 
consistent with our concern on the IR limit and the leading corrections. 
The equations of motion are 



an 



(A-l)r 5 



n" f" n' 2 
77 ~ 2f + 2rf 



n'f , 5f 2 /' 



nf 8f 2 2rf 



2(2A - l)r 



n 



n 



-(3C + 877) (f"f - \f 2 ^ - 1 [(5C + 14r?)(/ - l) 2 + 2(3( + 8r?)(/ - 1)] - £(/ - 1) 



a 5 
+ -r 2 
2 



and 



a(A- 1) 



2n" f" n'f f 



r> 



- -ir 7 + 4r 4 



n 



f nf'f 



n 



n 



TW ~ V- ) • < 10 » 



+ 2ar J - = -— -. 11 
/ n on 
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Here the equations of motion have been obtained by directly inserting the metric ansatz 
(J7|) into the action and then by variation of the action. Since the metric has good symme- 
tries, (JTQ]) - (|lip coincide with the Euler-Lagrange equations for static spherically symmetric 
objects, derived from the action (j3J) after assigning the projectability condition N = 1. 

Now we consider the vacuum solution for which £m = 0. Concerning the character of 
the solution, we consider the A = 1 case for tractability, which is also consistent with the 
recovery of GR in the IR limit. A newly proposed version of HL gravity, which is free 
from the unwanted scalar graviton, possesses an extra local U(l) symmetry to the foliation- 
preserving diffeomorphism naturally fixes A to be unity 0. With A = 1 the equation (fTTj) 
is reduced to /' = so that we write a constant solution 

/(r) = l + / >0. (12) 

Inserting ffl2|) into the equation ([TO]) , we obtain 



(rfy = -±J> + i /o + (5C + 14^ + 2(3C + 8,)/ 
2a a ar z 

The solution to this equation is 



o 2 M £ , M r s (5C + 14r ? )/ 2 + 2(3C + 8r / )/ 

n(r) = ±\ ~— r 2 + -f + , 14 

ba a r ar^ 



where r s is the integration constant which can be identified as r s = 2GM and M is the mass 
of black hole as in the Schwarzschild solution. Two comments are noted before analyzing 
the solution. First, if fo = 0, the solution reduces to that of GR even for £/a ^ 1 and in the 
presence of the higher spatial derivative terms. Second, the full set of field equations of the 
non-projectable version of HL gravity up to quartic spatial derivative terms were derived in 
j8|, Tl]]. The projectability condition N = N(t) is not compatible with those equations in 



general. Since the non-projectable version is a different theory from the projectable version 
of our interest, our analysis will be focused on the obtained solution ( 111]) . 

In the GR limit where £/a = 1 and a = £ = rj = 0, fo can be rescaled to be zero by an 
appropriate t-r mixing coordinate transformation. On the other hand, in HL gravity, such 
t-r mixing coordinate transformation is not allowed as a symmetry transformation. Thus, 
we have an additional integration constant /q. Let us examine the implications of this new 
integration constant. With a = 0, the solution is given by 



/(r) = l + / 0) n(r) = ±J£/ + — 4, ( 15 ) 

V a r r z 

where d = [(5( + 14r/)/ 2 + 2(3( + 8r])f ]/a. When £/a = 1, it looks the same as the 
Reissner-Nordstrom solution in the PG coordinates with the identification d = Gq 2 where q 
is the electric charge. When fo ^ 0, the l/r 2 -term acts like that of an electric charge. The 
important difference from the Reissner-Nordstrom solution of GR is that the coefficient of 
l/r 2 -term can have a positive value when £ or rj < 0. Though the quadratic curvature terms 
induce possibility of d = Gq 2 > 0, we will read the meaning of it in the context of GR. If 
we write down the Einstein equation, G^ v = —8irGT^, in terms of the metric (j7j) with (1121) . 
it becomes 



\rn 



- f - MGT\r 2 = f + 4vrG/o£ r V, (16) 
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where the last equality holds for the radial component of electrostatic field E r = F rt . Com- 
paring ( fl6|) with ( Tl3|) in the limit of a = r s = £/a — 1 = 0, we identify the energy density 
as 

Therefore, in order to obtain the same solution n(r) (115j) in GR, we need E% = q 2 /Anf^ 
which can be negative for d = Gq 2 < 0. In the context of GR, it is forbidden since this 
matter configuration violates the positive energy theorem. In conventional electromagnetism 
the effect of q 2 < may imply sign-flipped Coulomb force, attractive between same charges. 
However, in HL gravity, the solution ( T15]) with negative d is generic. It is obtained in the 
absence of matter, Cm — 0, but in the presence of quartic spatial derivative terms, the 
third term in the right-hand side of the equation f fl~3|) . Since violation of the positive energy 
theorem and modification of such Coulomb force are unphysical, detection of the signal of 
—d/r 2 > term suggests existence of the era of HL gravity. 

If fo cannot be scaled to zero in the metric of (I15p . another intriguing question is to 
address the effect of the constant term £/o/ a which is unphysical in GR and is not related 
to higher spatial derivatives. To see this explicitly, let us examine the geodesic equation of 
a test body. Let the constants of motion corresponding to the cyclic coordinates t and be 
E and £, respectively. Then radial geodesic equation can be written as 



^+|l + / i-^+*lll + ^)=B", (18) 



where the overdot denotes differentiation with respect to the proper time of the test body. 
We introduced the factor 5 = 1 — £/a measuring the deviation of the propagation speed 
from unity. From this equation, we obtain the orbit equation for u(<f)) = l/r(0) 

+ k5+-y = ^- 2 + ^u 2 + 2du\ (19) 

When we turn off the higher derivative terms d = and the mass r s = 0, the orbit equation 
reduces to a linear equation and the coefficient of n-term decides the allowed range of 
the angle 0. Solution of it is u(<p) oc cos(v / l + foS 0), and thus, non-zero fo leads to the 
surplus solid angle A = 27r(i/l + f Q 5 — 1) for f 5 > and the deficit solid angle A = 
27r(l — a/1 + fo5 ) for —1 < f Q 5 < 0. When f 5 = 0, the space has neither surplus nor deficit 
solid angles. 

In the context of GR (fl6l) . the constant metric solution, n = £fo/at, from (Tl3|) can be 
obtained by assuming the energy density, — T\ = —fo5/8irGr 2 . For the deficit solid angle 
with —1 < fn6 < 0, the energy density is positive and it corresponds to gravitating global 
monopole [LJ|. For the surplus solid angle with f 8 > 0, the energy density is negatively 
distributed everywhere. It is forbidden in GR since it violates the positive energy theorem. 
However, in HL gravity, this constant solution is a static solution attained in the absence 
of matter field, Cm = 0, and higher derivatives, ( = rj = 0. This result is also contrasted 
with the case of HL gravity with the detailed balance condition, in which a positive energy 
distribution of the electrostatic field of a point charge supports the geometry involving 
surplus solid angle [THI ] . 

Possible astrophysical effects of a deficit /surplus solid angle can easily be visualized in 
the case of 9 = ir/2, where an observer, a light source, and the the apex of deficit /surplus 
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solid angle are in the same plane. When a geometry with deficit solid angle is formed, the 
light from a star behind the apex propagates straight and arrives at a static observer who 



detects double images projected behind the source [16j. When a geometry with surplus solid 



angle is formed, a static observer tracking down the trajectory of the star experiences sudden 
disappearance of its image for the period proportional to the surplus angle and reappearance 
at a distant point over the apex of surplus angle [HI]. The present astronomical bound of 
angular resolution is about 200 micro- arcseconds (~ 10~ 9 radian) It means that, 

in principle, it can be observed if the graviton speed deviates away from unity satisfies 
\S\ = |1 — £/a| > 10 _9 7r//o- Though f is an undetermined free parameter in the present 
stage, the formula suggests 10~ 8 ~ 10~ 9 deviation of the graviton speed for f ~ 0(1), 
which looks extremely stringent for the RG flows for the ratio £/a in the vicinity of IR fixed 
point. 



IV. STATIC SPHERICALLY SYMMETRIC SOLUTION IN THE POINCARE 

COORDINATES 

In this section, we consider static spherically symmetric vacuum solution in non- 
projectable HL gravity for comparison. This is done by taking a static spherically symmetric 
metric in the Poincare coordinates 



ds 2 = -Nirfdi 2 + -^rrdr 2 + r 2 (d9 2 + sin 2 9d(f>) 2 . 
F{r) 

The metric ([7]) is connected to the metric (120]) by a coordinate transformation 



(20) 



I ^„ F dr. 



(21) 



with r, 9, unchanged and 



f = C 



N 2 '' 



n 2 = {C- N 2 ) 



N 2 '- 



(22) 



where C is a constant. Since this transformation is not a foliation-preserving diffeomorphism, 
the solutions found using the metric (120]) and the metric (J7|) will not be equivalent in HL 
gravity. 

The action (HI) written in terms of N and F is 



Sht, = 47r / dt dr 



r 2 N 



~^(rF' + F-l) + a 



(23) 



and the matter-field action with spherical symmetry is 

f r 2 N 
S M = 4vr J dt dr —^= C M (N, F, $). 



(24) 



8 



The equations of motion obtained from the variation of N and F are 



3C + 8 



2r 2 



log 



V F ,2 



2(C + 4r/) 



F'(F 



2(C + 2 V ) 



2 e 



tF' + F 



C 



1) + a 
dC M 



N 



dN 



N 



3C + 877 ^ 2(C + 4ry) ^ 
— 3 — * + 3 



r 



(3C + 877) 



F "--.( F - 1 

dC M N dC M 
dF F dN ' 



(25) 



(26) 



where the prime denotes the derivative with respect to r. Now we consider the vacuum 
solution of ( )25|) - (l26|) for which Cm = 0. The first equation involves F(r) only and the 
second equation determines N(r) from F(r) obtained from the first equation. Note that 
both equations are independent of a, (3, and 7 since the extrinsic curvature and the Cotton 
tensor vanish for the static spherically symmetric metric ( |20|) . When compared to (fTOj) and 
(fTTj) . a notable difference is that (125]) and (1261) are independent of A. 
We write the metric function F as 



where 



F(r) 



1 + r 2 [a — p(r)] 



k(C + 377) 



4(C + 3tj)' 

Then, from fl25l) . p(r) satisfies the equation 

(3C + 8r/)rV 2 + 4[4(C + 3r))p + |]r^ + 12[2(C + 3r/)p 2 
When both C and rj vanish, the well-known GR solution 

N 2 (r) = F(r) 



b\ = 0. 



1 - — r 2 - — 

6^ r 



(27) 
(28) 

(29) 
(30) 



is obtained, regardless of the coefficient of scalar curvature term, £. Here r s is the same inte- 
gration constant as in ( fl4l) . r s = 2GM. Note that when we have nonvanishing cosmological 
constant, a 7^ 0, the asymptotic behavior is determined by the coefficient of r 2 term. When 
higher derivative terms are introduced, its coefficient is modified and hence the asymptotic 
behavior becomes different from that of GR, irrespective of the assignment of the detailed 
balance condition. 

When the combination 3( + 877 vanishes, the equation (1291 is easily integrated and the 
solution is given by 



N 2 (r) = F(r) = 1 + ar 2 - 1 



2r_ 2 




(31) 



where £ = Let us consider the case a = 0. Then, we have a = and £ = £. For large r 
(r 3> (C r s) 1//3 ), we get an approximation 



(32) 
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which shows the usual behavior of the Schwarzschild black hole plus small corrections due 
to the higher derivative terms. On the other hand, at short distance, the position of event 
horizon is modified. When £ < 0, we have two event horizons at 



Tu — — 

H 2 



r e ± \/r? 



-0 



When C > 0, we have an event horizon at 

1 



= 2 I r s+ \/ r s +C 

For 3( + 8r] ^ 0, we solve (|29|) for p'(r) and obtain 

6 



(33) 



(34) 



P' 



where p = (p and 6 = 1 + 3r//£. When 6 = 0, this equation is easily integrated to give 



(35) 



J 2r 3 ' 



(36) 



For 6 7^ 0, we get 



/> 2 



1 + 2b (2 + p + 2 y/l - p - 2bp 2 ) ] (^M + 2y/l-p-2bp 2 



2-p + 2^1-p-2bp 2 



c 
j.6 



(37) 



where c is a constant. To get p(r) explicitly, we need to invert these equations. Assuming p 
is small, we can solve the equation perturbatively and obtain 



pas — 



4r 3 



Then up to this order 



F(r) 



1 I 2 ' s 

1 + ar 



Cr s 2 



N 2 



F{r) + 



4£r 4 
3aC(3C + 8r/K 2 



4^2 r 4 



(38) 

(39) 
(40) 



In general, the higher derivative terms give rise to the subleading corrections of order (r s /^r 3 . 
Thus, Eddington-Robertson parameters are same as those of GR, and it is hard to observe 
the macroscopic effect of them. However, higher spatial derivative terms modify the causal 
structure at short distance. 

In the previous and present sections, we obtain static vacuum solution of a projectable 
and a non-projectable version of HL gravity without the detailed balance condition at the 
IR regime, which are differentiated by two inequivalent metrics. For vanishing cosmological 
constant case, static spherically symmetric solution for the metric without nontrivial lapse 
function in the PG coordinates supports geometry of a surplus or deficit solid angle, however 
the other solution in the Poincare coordinates does not. Instead of such a long range effect, 
it involves a short distance correction like the change of event horizons. 
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V. CONCLUSION 

We obtained static spherically symmetric vacuum solutions of the recently proposed HL 
gravity, imposing the detailed balance condition only in the UV limit. Since HL gravity is 
intrinsically non-relativistic, the choice of a preferred frame is unavoidable. The choice of 
coordinate system is related to the choice of the frame. For static spherically symmetric 
vacuum solutions, we tried two coordinate systems, the PG coordinates and the Poincare 
coordinates. They are connected by the time-space mixing coordinate transformation which 
is not foliation-preserving. Thus, the solutions in two coordinate systems are distinct and 
physically inequivalent. The distinguishing feature of metrics in two coordinate systems is 
that the one in PG coordinates satisfies the projectability condition and the other does not. 

Even in the absence of higher derivative terms, the solution without matter distribution 
in the PG coordinates has an additional integration constant which leads to the geometry 
of a surplus or deficit solid angle. This long range effect can be constrained by smallness of 
the astronomical bound of angular resolution, and so do some parameters of HL gravity. In 
both solutions, the effect of higher spatial derivative terms is that they change the causal 
structures at short distances. But the differences from GR are not observationally significant 
in the macroscopic world for the reasonable choice of parameters. This is as expected because 
the higher spatial derivative terms are in general suppressed by the Planck scale or, if any, 
the new quantum gravity scale which may be smaller than the Planck scale but still out of 
our experimental reach. 

The solution in PG coordinate system may be stated as the static spherically symmetric 
vacuum solution of projectable HL gravity. For this solution we found that a new integration 
constant must be introduced in addition to mass, when £/a differs from unity or we have 
non- vanishing higher spatial derivative terms. This constant implies the existence of surplus 
or deficit angle for the former and the electric charge for the latter. The current astronomical 
bound on surplus/deficit angle can impose stringent constraint on the deviation of £/a from 
unity. One striking feature of the solution is that the electric charge can effectively be 
imaginary depending on the signature of the couplings of higher spatial derivative terms. 
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